Phase shift effects on the second moment and skewness of the field profiles obtained 

by the muon spin relaxation technique 
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Abstract. Recent high transverse field muon spin relaxation (TF-/iSR) experiments performed 
on Bi2212 single crystals show that the phase parameter extracted from individual histogram, an 
indication of the angle between the muon's initial polarization direction and the detector direction, is 
temperature dependent. This phase shift effect, probably due to the electronic or initial polarization 
instability at a time interval comparable with the muon spin precession period, will affect the second 
moment and skewness of the field profiles. The proper data analysis procedure is discussed to correct 
this phase shift effect, which is important on the quantitative analysis of the first order transition 
or melting transition in the Bi2212 mixed state. 
Keywords: Phase shift; Muon spin relaxation; Second moment; Skewness 
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I. INTRODUCTION 

In a /iSR experiment ||l| the spin-polarized muons are 
assumed to stop at random positions within the sam- 
ple and precess along the local (internal) magnetic fields 
B{r). The procession frequency is 7^ -B, where 7^ = 135.5 
MHz/T is the gyromagnetic ratio of the muon. If the 
local fields which the muons feel are not unique i.e. there 
is a distribution of magnetic field f{B), consequently 
the muons will precess with different frequencies, which 
causes a depolarization of the muons ensemble's spins. 
The second moment of the field distribution f{B) indi- 
cates how fast the muons depolarize thus is called the 
relaxation rate {a or A). There are several ways to ex- 
tract the relaxation rate: we can assume the proper form 
of the depolarization function P{t) to fit the /iSR spectra 
(histogram or asymmetry) in the time domain to extract 
the relaxation rate parameter in the fitting procedure; or 
we can calculate the second moment of the field distribu- 
tion f{B) to get the relaxation rate a since the Cosine 
Fourier Transformation (CFT) of /iSR spectra indicates 
the magnetic field distribution inside the sample. pU2] 
The second moment calculated from the discretized field 
distribution data depends slightly on the field channel 
width A, there exists the following relations: p] 
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where a is the true relaxation rate and Udata is the re- 
laxation rate calculated from the discretized field profile. 
Simulation shows that as long as the channel width A is 
much smaller than the second moment of the field profile, 
(^data is a good estimate of a. H For a ^SR spectra in 



a 10 /IS time window, the ideal minimum relaxation rate 
we can extract from data is 0.628 ^.s~^ . If we consider 
the channel width effect mentioned above, the resolution 
of the relaxation rate will still be larger. 
The local fields f{B) can either be intrinsic, as they are 
for ordered magnets and spin glasses, or induced by an 
external field, as for the vortex lattice of a type II super- 
conductor formed in the external magnetic fields. Usually 
the distribution f{B) is dependent on the flux distribu- 
tion associated with a single vortex line, as well as the 
arrangement and dynamics of the vortex lines. Note for 
single vortex, the magnetic field has a component normal 
to the vortex axis, however, the average transverse com- 
ponent vanishes due to the large amount of contributions 
from the entire FLL H. For general calculation of the 
magnetic field distribution, see the monograph of Greer 
and Kossler |^. 



II. PHASE-SHIFTED /iSR SPECTRA 

In most cases, there is no difficulty to obtain the true 
field profiles through CFT from low transverse field /iSR 
data. However, when the applied magnetic field is so high 
(several Tesla) that one precession period is only several 
times of the channel width, we need to study the effect of 
time-shifted or phase-shifted P{t) (due to the instabihties 
of electronics or initial polarization) on the interpretation 
of f(B). Assume P(t) = Gt(0 cos(woi) where Grit) is 
the non-oscillating depolarization function part and can 
be exponential, stretched exponential, Gaussian or other 
types. Use the notations 
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We can easily derive the Fourier Transform (FT) of the 
time-shifted (or phase-shifted) P{t) as follows: 



Pc{uj) = / GT{t)cos{ujt) dt 
Ps{uj) = / GT{t)sm{ujt) dt 



Psc{i^) = I Grit + St) cos[ijjo{t + St)] cos{ujt) dt 



cosjuJoSt) sinjuJoSt) 

Pc(UJo - W) Ps[UJo - UJ) (4) 



underlying changes of the vortex structure (the one-to- 
one mapping still needs to be studied further) , is defined 
from the 3rd and 2nd moments of the field line shape as 
follows: 
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smjupSt) cosjupSt) 
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Note in the expression of Psc(w), the 1st term is sym- 
metric about uq, however, there is a sign change in the 
2nd term when uj goes across loq- This sign change was 
seen in our Bi2212 high transverse field /xSR data |^ at 
various temperatures. 

To construct the true field profile PdoJo — <^) from the 
experimental data Psd^) and Pssijjj), we can do the fol- 
lowing transformation, 

Psci^) cos{ujoSt) - Pss(w) sm{ujoSt) = -Pc{loo - w) (6) 

The result of this phase correction on one of the field 
profiles in Bi2212 is shown in figure 1 where we can see 
the proper field profile is obtained. 
The measurement of the internal magnetic field distribu- 




2297 2299 2301 

0) (MHz) 



FIG. 1. The field profile on the ab basal plane of Bi2212 
obtained before (the dashed line) and after (the solid line) the 
phase correction as discussed in the text. 

tion in the vortex of type II HTSC allows an investigation 
of the details of the static or dynamic flux structure [g , 
such as the effects of random pinning g and the dimen- 
sionality/melting of the vortex ^!^. To quantify and 
characterize the shape of the field distribution (field pro- 
file), skewness factor a, the variation of which reflects the 
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where M„ = J B"f{B)dB. and is the nth moment of 

fiB). 

Normally we obtain the field profile cither from asymme- 
try plot or from individual histogram then calculate the 
second moment and skewness of the field profile to probe 
the possible phase transitions of the vortex matters. In 
the analysis of high transverse field /iSR data on Bi2212, 
we often find that the phase parameter in the depolariza- 
tion function changes with temperature and thus needs 
to be corrected precisely. 



III. CALCULATION RESULTS 

To evaluate the effects of the phase shift (p on the sec- 
ond moment and skewness calculations of the field profile, 
we generate a depolarization function (asymmetry plot) 

P(t) = e 2— cos(a;i -I- 4>) where a — 2 fj.s~ , w = 20 
Mrad/s (w can be assumed much higher, the final results 
are similar). These parameters are chosen so that there 
are no evident visual changes in the field profiles when 
the phase shift varies from to 0.5 (^^30"). The cal- 
culated (j) dependent second moment (ideally the second 
moment should be 2 ^J,s^^) and skewness plot is shown 
in the figure 2. 

We can see clearly that the second moment and skew- 
ness are very sensitive on the phase shift (f). A change of 
initial phase from to 0.5 (~30°) causes about 50 percent 
change on the relaxation rate and an increase of skewness 
from to 2. This result requests us to extract the phase 
parameter accurately after fitting the /iSR asymmetry 
plot. To correct this phase-shift effect, we can follow the 
modified CFT procedure as discussed in the previous sec- 
tion to obtain the true field profiles, which is also the case 
when we analyze the individual histogram. 
We now study the asymmetry A(t) which is used to ob- 
tain the field profile. We know A{t) is usually composed 
by two histograms coming from two oppositely (tt out 
of phase) located positron detectors. However, in real- 
ity, especially for high precession frequency data, the two 
histograms are usually found not exactly tt out of phase. 
To consider this phase deviation (j) from tt, we can do 
the following analytical analysis. As usual, we assume 
the phase-shifted histograms are (we assume Gaussian 
depolarization here, we can also use other depolarzation 
function instead) 



IV. CONCLUSION 
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FIG. 2. The phase shift effects on the second moment and 
skewness of the field profile, detailed procedure is discussed 
in the text. 



Nu{t) = Noe--[l + Aoe-'^^ cos(wt)] (9) 

Noit) = A^oe"- [1 - Aoe-'^ cos(wt + 0)] (10) 

The asymmetry A{t) can be expressed as 
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From equation 11, wc can sec that, to the first order, 
the relaxation function is composed of two depolarized 
harmonics with frequencies uj and 2uj. The asymmetry 
amplitude for the uj harmonic is j4oCOs(^), which is re- 
duced by a factor of cos(^) from the original amplitude 
Aq , there is a phase shift ^ , which causes the changes on 
the second moment and skewness of the field profile as we 

have discussed before; since Aq is usually less than 0.4, 

^2 
the asymmetry amplitude for the 2uj harmonic -^ sin (p 

is much less than Aq cos(^), the amplitude of the oj har- 
monic. All these conclusions have been verified through 
further simulation. 



In order to obtain the proper field profiles through the 
Cosine Fourier Transform (CFT) on the /^SR data and 
characterize them quantitatively and correctly, wc need 
to extract the initial phase parameter precisely from indi- 
vidual histogram or asymmetry plot and correct the pos- 
sible phase shift effect through a proposed data analysis 
procedure in order to obtain the correct second moment 
(or relaxation rate) and skewness of the field profile. 
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